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We propose to realize a mass gap in QCD by not imposing the transversality condition on the 
full gluon self-energy, while preserving the color gauge invariance condition for the full gluon prop- 
agator. This is justified by the nonlinear and nonperturbative dynamics of QCD. None of physical 
observables/processes in low-energy QCD will be directly affected by such a temporary violation of 
color gauge invariance/symmetry. No truncations/approximations and no special gauge choice are 
made for the regularized skeleton loop integrals, contributing to the full gluon self-energy, which 
enters the Schwinger-Dyson equation for the full gluon propagator. In order to make the existence 
\ j^ ' of a mass gap perfectly clear the corresponding subtraction procedure is introduced. All this allows 

f^ I one to establish the general structure of the full gluon propagator and the corresponding gluon 

C^ . Schwinger-Dyson equation in the presence of a mass gap. It is mainly generated by the nonlinear 

^S| ' interaction of massless gluon modes. The physical meaning of the mass gap is to be responsible for 

the large-scale (low-energy /momentum), i.e., nonperturbative structure of the true QCD vacuum. 
In the presence of a mass gap two different types of solutions for the full gluon propagator are 
I ■ possible. The massive solution leads to an effective gluon mass, which depends on the gauge-fixing 

^— ( ' parameter explicitly. This solution becomes smooth at small gluon momentum in the Landau gauge. 

w-> , The general iteration solution is always severely singular at small gluon momentum, i.e., the gluons 

remain massless, and this does not depend on the gauge choice. 

,-j. ' PACS numbers: 11.15.Tk, 12.38.Lg 
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Tjj- . I. INTRODUCTION 

o 

^y • Today there is no doubt left that color confinement and other dynamical effects, such as spontaneous breakdown 

^r ^ of chiral symmetry, bound-state problems, etc., being essentially nonperturbative (NP) effects, are closely related to 

,!-{ . the large-scale (low-energy/momentum) structure of the true QCD ground state and vice- versa [l|, Q (and references 

I ' therein). The perturbation theory (PT) methods in general fail to investigate them. If QCD itself is a confining 

C^ theory then a characteristic scale has to exist. It should be directly responsible for the above-mentioned structure of 

~5 ' the true QCD vacuum in the same way as Aqcd is responsible for the nontrivial perturbative dynamics there (scale 

violation, asymptotic freedom (AF) Q). 

The Lagrangian of QCD [3, |j| does not contain explicitly any of the mass scale parameters which could have a 
physical meaning even after the corresponding renormalization program is performed. The main goal of this paper is 
to show how a characteristic scale (the mass gap, for simplicity) responsible for the NP dynamics in the infrared (IR) 
C^ ' region may explicitly appear in QCD. This becomes an imperative especially after Jaffe and Witten have formulated 
their theorem "Yang-Mills Existence And Mass Gap" [3]. We will show that the mass gap is dynamically generated 
mainly due to the nonlinear (NL) interaction of massless gluon modes. 

The propagation of gluons is one of the main dynamical effects in the true QCD vacuum. It is described by the 
corresponding quantum equation of motion, the so-called Schwinger-Dyson (SD) equation [ll] (and references therein) 
for the full gluon propagator. The importance of this equation is due to the fact that its solutions reflect the quantum- 
dynamical structure of the true QCD ground state. The color gauge structure of this equation is the main subject of 
our investigation in order to find a way how to realize a mass gap in QCD. Also we will discuss at least two possible 
types of solutions of the gluon SD equation in the presence of a mass gap, making no approximations/truncations 
and no special gauge choice for the skeleton loop integrals contributing to it. So they can be considered as the 
generalizations of the explicit solutions because the latter ones are necessarily based on the above-mentioned specific 
approximations/truncations schemes. 
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II. QED 

It is instructive to begin with a brief explanation why a mass gap does not occur in quantum electrodynamics 
(QED). The photon SD equation can be symbolically written down as follows: 

D{q) = D''{q) + D^{q)Il{q)D{q), (2.1) 

where we omit, for convenience, the dependence on the Dirac indices, and D^{q) is the free photon propagator. Il{q) 
describes the electron skeleton loop contribution to the photon self-energy (the so-called vacuum polarization tensor) . 
Analytically it looks 

n(g) EE n^.(g) = -g^ J ^Tr[j^S{p - g)r.(p - g, q)S{p)], (2.2) 

where S{p) and T^{p — q,q) represent the full electron propagator and the full electron-photon vertex, respectively. 
Here and everywhere below the signature is Euclidean, since it implies gi — > when g^ — > and vice-versa. This 
tensor has the dimensions of a mass squared, and therefore it is quadratically divergent. To make the formal existence 
of a mass gap (the quadratically divergent constant, so having the dimension of a mass squared) perfectly clear, let 
us now, for simplicity, subtract its value at zero. One obtains 

ff (q) = xi^M = np.(g) - n^.(o) = n^.(g) - s^,a^{x). (2.3) 

The explicit dependence on the dimensionless ultraviolet (UV) regulating parameter A has been introduced into the 
mass gap A^(A), given by the integral (2.2) at q^ — 0, in order to assign a mathematical meaning to it. In this 
connection, a few remarks are in order in advance. The dependence on A (when it is not shown explicitly) is assumed 
in all divergent integrals here and below in the case of the gluon self-energy as well (see next section). This means 
that all the expressions are regularized (including photon/gluon propagator), and we can operate with them as with 
finite quantities. A should be removed on the final stage only after performing the corresponding renormalization 
program (which is beyond the scope of the present investigation, of course). Whether the regulating parameter A has 
been introduced in a gauge- invariant way (though this always can be achieved) or not, and how it should be removed 
is not important for the problem if a mass gap can be " released/liberated" from the corresponding vacuum. We will 
show in the most general way (not using the PT and no special gauge choice will be made) that this impossible in 
QED and might be possible in QCD. 

The tensor structure of the subtracted photon self-energy can be written as follows: 

K-^iq) = T^.(<?)g'n^(<z') + qt^qAqMiq^), (2.4) 

where both invariant functions 11^ (q^) at n = 1,2 are, by definition, dimensionless and regular at small q^, since 
n*(0) = 0; otherwise they remain arbitrary. From this relation it follows that H'^iq) = 0{q^), i.e., it is always of the 
order q^ . Also, here and everywhere below 

T^i^iq) = S^iy - qf^q,y/q^ = Sf^,. - Lf^„{q). (2.5) 

Taking into account the subtraction (2.3), the photon SD equation becomes 

D{q) - D°{q) + D"{q)W{q)D{q) + D" (q) A^X) D (q) . (2.6) 

Its subtracted part can be summed up into the geometric series, so one has 

D{q) = D"{q) + D"{q)A'{X)D{q), (2.7) 

where the modified photon propagator is 



D\q) - i_ns(g^j^O(g) - ^°(9) + D'\qW{q)D\q) D%qW {q)D\qW {q)D^\q) + ... . (2.8) 

Since Il-^{q) = 0{q^) and D^{q) ^ (<z^)~^, the IR singularity of the modified photon propagator is determined by the 
IR singularity of the free photon propagator, i.e., D^{q) = 0{D^{q)) with respect to the behavior at small photon 
momentum (in Eqs. (2.6), (2.7) and (2.8) we again omit the tensor indices, for simpHcity). 

Similar to the subtracted photon self-energy, the photon self-energy (2.2) in terms of independent tensor structures 
is 

n,,.(<z) - T^.q2ni(q2) + q^g,n2(g2), (2.9) 

where both invariant functions n„(q^) at n= 1,2 are dimensionless and remain arbitrary. Due to the transversality 
of the photon self-energy 

q^J^^.u{q)^qu^^.u{q) = Q. (2.10) 

which comes from the current conservation condition in QED, one then has Ii2{q^) — 0, i.e., this tensor should be 
purely transversal. On the other hand, from the subtraction (2.3) and transversality condition (2.10) it follows that 

n^(g2) = -(A2(A)/g2). (2.11) 

However, this is impossible, since 112(9^) is a regular function of q^ , so the mass gap should be zero and consequently 
n|(q^) = as well, i.e., 

n^(g2) = A2(A) = 0. (2.12) 

This means that the photon self-energy coincides with its subtracted counterpart and both of them are purely transver- 
sal. Moreover, this means that the photon self-energy does not have a pole in its invariant function ni(g^) = Ilf (g^). 
As mentioned above, in obtaining these results neither the PT has been used nor a special gauge has been chosen. 
So there is no place for quadratically divergent constants in QED, while logarithmic divergence still can be present 
in the invariant function ni(g^) = nf(q^). It is to be included into the electric charge through the corresponding 
renormalization program (for these detailed gauge-invariant derivations explicitly done in lower order of the PT see 
Refs. 0,i,0,ii). 

In fact, the current conservation condition (2.10) lowers the quadratical divergence of the corresponding integral 
(2.2) to a logarithmic one. That is the reason why in QED logarithmic divergences survive only. Thus in QED 
there is no mass gap and the relevant photon SD equation is shown in Eq. (2.8), simply identifying the full photon 
propagator with its modified counterpart. In QED we should replace n(g) by its subtracted counterpart n'*((7) 
from the very beginning (n(g) — » n*(q)), totally discarding the quadratically divergent constant A^(A) from all the 
equations and relations. The current conservation condition for the photon self-energy (2.10), i.e., its transversality, 
and for the full photon propagator q^q^D^^i^q) = i^, where ^ is the gauge-fixing parameter, are consequences of 
gauge invariance. They should be maintained at every stage of the calculations, since the photon is a physical 
state. In other words, at all stages the current conservation plays a crucial role in extracting physical information 
from the S'-matrix elements in QED. For example, if some QED process includes the full photon propagator, then 
the corresponding 5-matrix element is proportional to the combination j'^{q)D^^{q)J2{q). The current conservation 
condition jf (g)?^ = J2iq)q'^ ~ ^ implies that the unphysical (longitudinal) component of the full photon propagator 
does not change the physics of QED, i.e., only its physical (transversal) component is important. In its turn this means 
that the transversality condition imposed on the photon self-energy is important, since 11^^ (g) itself is a correction to 
the amplitude of the physical process, for example such as electron-electron scattering. 

III. QCD 

Due do color confinement in QCD the gluon is not a physical state. Still, color gauge invariance should also be 
preserved, so the color current conservation takes place in QCD as well. However, in this theory it plays no role 



in the extraction of physical information from the S'-matrix elements for the corresponding physical processes and 
quantities. So in QCD there is no such physical amplitude to which the gluon self-energy may directly contribute (for 
example, quark-quark/ antiquark scattering is not a physical process). The lesson which comes from QED is that if 
one preserves the transversality of the photon self-energy at every stage, then there is no mass gap. Thus, in order to 
realize a mass gap in QCD, our proposal is not to impose the transversality condition on the gluon self-energy, but 
preserving the color gauge invariance condition for the full gluon propagator (see below). As mentioned above, no 
QCD physics will be directly affected by this. So color gauge symmetry will be violated at the initial stage (at the 
level of the gluon self-energy) and will be restored at the final stage (at the level of the full gluon propagator). 

A. Gluon SD equation 

The gluon SD equation symbolically is 

D^.u{q) = DIM + DlM^n^M, D)D,,{q), (3.1) 

where D'^^{q) is the free gluon propagator. IIpcy{q;D) is the gluon self-energy, and in general it depends on the full 
gluon propagator due to the non-Abelian character of QCD (see below). Thus the gluon SD equation is highly NL, 
while the photon SD equation (2.1) is a linear one. In what follows we omit the color group indices, since for the gluon 
propagator (and hence for its self-energy) they are reduced to the trivial (5-function, for example D'^'Kq) = D^^{q)S°-^. 
Also, for convenience, we introduce i into the gluon SD equation (3.1). 
The gluon self-energy Ilp^{q; D) is the sum of a few terms, namely 

n,,(g; D) - -UIM - KM) + KoiD) + n(i),,(g; D) + n(2),,(g; D) + n'(2),,(g; D), (3.2) 

where 11^^. [q) describes the skeleton loop contribution due to quark degrees of freedom (it is an analog of the vacuum 
polarization tensor in QED, see Eq. (2.2)), while Ii^I^{q) describes the skeleton loop contribution due to ghost 
degrees of freedom. Both skeleton loop integrals do not depend on the full gluon propagator D, so they represent 
the linear contribution to the gluon self-energy H* (13) represents the so-called constant skeleton tadpole term. 
n(i)po-((7; D) represents the skeleton loop contribution, which contains the triple gluon vertices only. Ti^2)pa{(l'i D) and 
n'/2-, (g; D) describe topologically independent skeleton two- loop contributions, which combine the triple and quartic 
gluon vertices. The last four terms explicitly contain the full gluon propagators in different powers, that is why they 
form the NL part of the gluon self-energy. The explicit expressions for the corresponding skeleton loop integrals |10| 
(in which the corresponding symmetry coefficients can be included) are of no importance here. Let us note that like 
in QED these skeleton loop integrals are in general quadratically divergent, and therefore they should be regularized 
(see remarks above and below). 

B. A temporary violation of color gauge invariance/symmetry (TVCGI/S) 

The color gauge invariance condition for the gluon self-energy (3.2) can be reduced to the three independent 
transversality conditions imposed on it. It is well known that the quark contribution can be made transversal 
independently of the pure gluon contributions within any regularization scheme which preserves gauge invariance, for 
example such as the dimensional regularization method (DRM) [3, [j, S l3) llll • So, one has 

(?,n^„(g) = g.n^,(g)=0, (3.3) 

indeed. In the same way the sum of the gluon contributions can be done transversal by taking into account the ghost 
contribution, so again one has 
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Il(i)pMD)+ n(2)p, (g; i?) + n'(2)p, (g; D) - Iif„ {q) =0. (3.4) 



The role of ghost degrees of freedom is to cancel the unphysical (longitudinal) component of gauge bosons in every 
order of the PT, i.e., going beyond the PT and thus being general. The previous relation just demonstrate this, since 
it contains the corresponding skeleton loop integrals. 



However, there is no such regularization scheme (preserving or not gauge invariance) in which the transversality 
condition for the constant skeleton tadpole term could be satisfied, i.e., qpIl*p^{D) = qpSp„A'^{D) = q^A'^{D) ^ 0, 
indeed. This means that in any NP approach the transversality condition imposed on the gluon self-energy may not 
be valid, i.e., in general 

qpHp^iq; D) - q,Wp„{q; D) / 0. (3.5) 

In the PT, when the full gluon propagator is always approximated by the free one, the constant tadpole term is set to 
be zero within the DRM ^Ei JJj]j i-e., H* (1)°) — 0. So in the PT the transversality condition for the gluon self-energy 
is always satisfied. 

The relation (3.5) justifies our proposal not to impose the transversality condition on the gluon self-energy. The 
special role of the constant skeleton tadpole term in the NP QCD dynamics should be emphasized. It explicitly 
violates the transversality condition for the gluon self-energy (3.5). The second important observation is that now 
ghosts themselves cannot automatically provide the transversality of the gluon propagator in NP QCD. However, this 
does not mean that we need no ghosts at all. Of course, we need them in other sectors of QCD, for example in the 
quark-gluon Ward-Takahashi identity, which contains the so-called ghost-quark scattering kernel explicitly [3| . 

C. Subtractions 

As we already know from QED, the regularization of the gluon self-energy can be started from the subtraction its 
value at the zero point (see, however, remarks below). Thus, quite similarly to the subtraction (2.3), one obtains 

n%{q; D) = n,,(g; D) - Hp,(0; D) - Hp,(g; D) - 5p,A\\- D). (3.6) 

Let us remind once more that for our purpose, namely to demonstrate a possible existence of a mass gap A^(A;£') 
in QCD, it is not important how A has been introduced and how it should be removed at the final stage. The mass 
gap itself is mainly generated by the nonlinear interaction of massless gluon modes, slightly corrected by the linear 
contributions coming from the quark and ghost degrees of freedom, namely 

A2(A; D) = n\D) + ^ H"(0; D) = A^D) + ^ Al{0; D), (3.7) 

a a 

where index "a" runs as follows: a — —q, —gh, 1, 2, 2', and — q, —gh mean that both terms enter the above-mentioned 
sum with minus sign (here, obviously, the tensor indices are omitted). In these relation all the divergent constants 
n*(D) and n''(0; D), having the dimensions of a mass squared, are given by the corresponding skeleton loop integrals 
at q^ = 0. Thus these constants summed up into the mass gap squared (3.8) cannot be discarded like in QED, since 
the transversality condition for the gluon self-energy is not satisfied, see Eq. (3.5). In other words, in QCD in general 
the quadratical divergences of the corresponding loop integrals cannot be lowered to logarithmic ones, and therefore 
the mass gap (3.7) should be explicitly taken into account in this theory. The transversality condition for the gluon 
self-energy can be satisfied partially, i.e., if one imposes it on quark and gluon (along with ghost) degrees of freedom 
as it follows from above. Then the mass gap is to be reduced to n*(Z?), since all other constants Il°'{0;D) can be 
discarded in this case (see Eq. (3.7)). However, we will stick to our proposal not to impose the transversality condition 
on the gluon self-energy, and thus to deal with the mass gap on account of all possible contributions. 
The subtracted gluon self-energy 

U;Aq;D)^U%q;D)=J2u:{q;D) (3.8) 

a 

is free from the tadpole contribution, because Ilf{D) — Hf(I?) — Ilt{D) = 0, by definition, at any D, while in the 
gluon self-energy it is explicitly present through the mass gap (see Eqs. (3.7) and (3.6)). The general decomposition 
of the subtracted gluon self-energy into the independent tensor structures can be written down as follows: 

Ii;,iq; D) = Tp„{q)q^n{q^-D) + qpqMi^; D), (3.9) 



where both invariant functions n((7^; D) and n(g^; D) are dimensionless and regular at small q^ . Since the subtracted 
gluon self-energy does not contain the tadpole contribution, we can now impose the color current conservation condition 
on it, i.e., to put 

gpn^,(<z;i?)-g.n^,(g;I?)-0, (3.10) 

which implies n(<7^; D) = 0, so that the subtracted gluon self-energy finally becomes purely transversal 

n%{q;D)=Tp„{q)q^n{q^;D), (3.11) 

and it is always of the order q^ at any D, since the invariant function I\-{q^; D) is regular at small q^ at any D. Thus 
the subtracted quantities are free from the quadratic divergences, but logarithmic ones can be still present in Tl{q^] D) 
like in QED. 

D. General structure of the gluon SD equation 

Our strategy is not to impose the transversality condition on the gluon self-energy in order to realize a mass gap 
despite whether or not the tadpole term is explicitly present. To show that this works, it is instructive to substitute 
the subtracted gluon self-energy (3.9) (and not its transversal part (3.11)) into the initial gluon SD equation (3.1), 
on account of the subtraction (3.6). Then one obtains 

D^.u{q) - D%{q) + Dlp{q)i[Tp,{q)q^n{q^;D) + q^q^^^; D)]DUq) + ^".('Z)*A2(A; D)D,,{q). (3.12) 

Let us now introduce the general tensor decompositions of the full and auxiliary free gluon propagators Df^u{q) — 
l[Tf,,{q)d{q^) + Lp,{q)di{q^)]{l/q^) and 

D%{q) - ^K^iq) + L^,{q)do{q^)]{l/q^), (3.13) 

respectively. The form factor do{q'^) introduced into the unphysical part of the auxiliary free gluon propagator £)° (q) 
is needed in order to explicitly show that the longitudinal part of the subtracted gluon self-energy Tl{q^; D) plays no 
role. The color gauge invariance condition imposed on the full gluon propagator 

q,,quDf,i,{q) =i^, (3.14) 

implies di{q'^) ~ ^, so that the full gluon propagator becomes 

Dp^q) = I {T^,{q)d{q^) + CL^Aq)} \- (3.15) 

Substituting all these decompositions into the gluon SD equation (3.12), one obtains 

'^''^ = i + n(,-z?) + (A^(A;^)/.^)' ^^-''^ 

and 

daiq^) = -~ ^ . (3.17) 

l-e[n((z2;D) + (A2(A;i?)/q2)] 

However, the auxiliary free gluon propagator defined in Eqs. (3.13) and (3.17) is to be equivalently replaced as follows: 



A2(A;D) 



DIM) =^ KM + ^iL^MdM) ^{q';D) + \ ' ^, (3.18) 
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where D^^{q) in the right-hand-side is the standard free gluon propagator, i.e., 

D%{q) = « {T^Aq) + ^L,,{q)} ^- (3.19) 

Then the gluon SD equation in the presence of the mass gap (3.12), on account of the exphcit expression for the 
auxihary free gluon form factor (3.17), and doing some tedious algebra, is also to be equivalently replaced as follows: 

+ DlM^^^'i^^-.mD.M+^eL^M^^^- (3.20) 

Here and below D^ {q) is the free gluon propagator(3.19). The gluon SD equation (3.20) does not depend on do{q'^) 
and n((j^; D), i.e., they played their role and then retired from the scene. So, our derivation explicitly shows that the 
longitudinal part of the subtracted gluon self-energy Il{q'^;D) plays no role and can be put to zero without loosing 
generality, and thus making the subtracted gluon self-energy purely transversal in accordance with Eq. (3.11). 
Using now the explicit expression for the free gluon propagator (3.19) this equation can be further simplified to 

D^..{q) - D%{q)-Tf,Aq)[^{q^;D) + ^lJ^^'^dUq)- (3.21) 

It is easy to check that the full gluon propagator satisfies the color gauge invariance condition (3.14), indeed. So the 
full gluon propagator is the expression (3.15) with the full gluon form factor given in Eq. (3.16), which obviously 
satisfies Eq. (3.21). The only price we have paid by violating color gauge invariance is the gluon self-energy, while 
the full and free gluon propagators and the subtracted gluon self-energy always satisfy it. Let us emphasize that the 
expression for the full gluon form factor shown in the relation (3.16) cannot be considered as the formal solution for the 
full gluon propagator, since both the mass gap A^(A; D) and the invariant function Il{q^; D) depend on D themselves. 
Here it is worth noting in advance that from above it is almost clear that if one begins with the UV renormalization 
program, then the information on the mass gap will be totally lost. In this case instead of the regularized gluon 
self-energy its subtracted regularized counterpart comes into the play. In other words, in the PT limit A^(A; _D) = 
one recovers the standard gluon SD equation, and the gluon self-energy coincides with its subtracted counterpart like 
in QED. 

Thus, we have established the general structure of the full gluon propagator (see Eqs. (3.15) and (3.16)) and the 
corresponding gluon SD equation (3.21) (which is equivalent to Eq. (3.12)) in the presence of a mass gap. 

IV. MASSIVE SOLUTION 

An immediate consequence of the explicit presence of the mass gap in the full gluon propagator is that a massive- 
type solution for it becomes possible. In other words, in this case the gluon may indeed acquire an effective mass. 
From Eq. (3.16) it follows that 

^'^^'^'^ = '^Tq^i^^o + mMY ^^'^^ 

where instead of the dependence on D the dependence on ^ is explicitly shown. The full gluon propagator (3.15) 
may have a pole-type solution at the finite point if and only if the denominator in Eq. (4.1) has a zero at this point 
q^ = —m? (Euclidean signature), i.e., 

-m2-m2H(-m2;O + A2(A,O = 0, (4.2) 

where mi = m^{X,^) is an effective gluon mass, and the previous equation is a transcendental equation for its 
determination. Excluding the mass gap, one obtains that the denominator in the full gluon propagator becomes 
q"^ + (j^H((j^; ^) -I- A^(A, C) = 9^ + rnl + q^^iq^',0 + "^gn(— m^; ^). Let us now expand Il{q'^-,S,) in a Taylor series near 
ml: 



n(g2; e) = n(-ml; + {q^ + m2)n'(-m2; + o[{q' + m^)^) . (4.3) 

Substituting this expansion into the previous relation and after doing some tedious algebra, one obtains q^ + m'i + 
g2n(g2;^) + m2n(-m2;e) = (<z2 + m2)[l + n(-m2;e)-m2n'(-m2;e)][l + n«(g2;^)], where n«(g2;^)=0 at (72 = _,Ti2 
(otherwise it remains arbitrary). Thus the full gluon propagator (3.15) now looks 

^-^('^) =*^-''('')(<z2+mp[l+V(g2;m2)] + ^^^^^^ ('^) ^ ' (4.4) 

where, for future purpose, in the invariant function II^{q^; mi) instead of ^ we introduced the dependence on the gluon 
effective mass squared m^ which depends on ^ itself. The gluon renormalization constant is Z^ = [I + n(— m^; ^) — 
m'^Il'{~m^;^)]~^. In the formal PT limit A^(A,^) = 0, an effective gluon mass is also zero, TOg(A,^) = 0, as it 
follows from Eq. (4.2). So an effective gluon mass is the NP effect. At the same time, it cannot be interpreted as the 
"physical" gluon mass, since it remains explicitly gauge-dependent quantity. The gluon renormalization constant in 
this limit becomes a standard one, namely [1 + 11(0; ■C)]"^- The massive-type solution (4.4) becomes smooth in the IR 
{q^ -^ 0) in the Landau gauge ^ = only (the ghosts now cannot guarantee the cancellation of the longitudinal part 
of the full gluon propagator as mentioned above) . In this connection let us point out that Landau gauge smooth (even 
vanishing in the IR) gluon propagator at the expense of more singular (than the free one) in the IR ghost propagator 
has been obtained and discussed (see, for example Refs. [l^jllJI ^^^ references therein). As mentioned above, however, 
these results are necessarily based on different approximations/truncations for the skeleton loop integrals contributing 
to the gluon self-energy. 

V. ITERATION SOLUTION 

In order to perform a formal iteration of the gluon SD equation (3.21), much more convenient to address to its 
"solution" for the full gluon form factor (3.16), nevertheless, and rewrite it as follows: 

d{q^) = 1 - \u(/;d) + ^!(Ml d{q^) = 1 _ p(q^;d)d{q'), (5.1) 

L q^ i 

i.e., in the form of the corresponding transcendental (i.e., not algebraic) equation suitable for the formal nonlinear 
iteration procedure. Here we replace the dependence on D by the equivalent dependence on d. For future purposes, 
it is convenient to introduce short-hand notations as follows: 

n(g2; d = d(0) + d(i) + d(2) + ... + dM + ...) = n„,(g2) ^ [P,„(g2) „ (AlJq% (5.2) 

In these relations A^ are the auxiliary mass squared parameters, while A^ is the mass gap itself (see, however, remarks 
in Conclusions) . The dimensionless constants Cm via the corresponding subscripts depend on which iteration for the 
gluon form factor d is actually done. They may depend on the dimensionless coupling constant squared g^, as well 
as on the gauge-fixing parameter ^. We also introduce the explicit dependence on the dimensionless finite (slightly 
different from zero) subtraction point a, since the initial subtraction at the zero point may be dangerous [3J. The 
dependence of A^ on all these parameters is not shown explicitly, and if necessary can be restored any time. Let us 
also remind that all the invariant functions Ilm{q^) are regular at small q^. If it were possible to express the full gluon 
form factor d{q'^) in terms of these quantities then it would be the formal solution for the full gluon propagator. In 
fact, this is nothing but the skeleton loops expansion, since the regularized skeleton loop integrals, contributing to the 
gluon self-energy, have to be iterated. This is the so-called general iteration solution. No truncations/approximations 
and no special gauge choice have been made. This formal expansion is not a PT series. The magnitude of the coupling 
constant squared and the dependence of the regularized skeleton loop integrals on it is completely arbitrary. 

It is instructive to describe the general iteration procedure in some details. Evidently, d^°^ = 1 and doing the first 
iteration in Eq. (5.1), one thus obtains 

d{q^) = 1 - Poiq') + ... = 1 + d«(g2) + ..., (5.3) 



where obviously ^(^^(q^) = —Po{q'^). Doing the second iteration, one obtains 

d{q^) = 1 - Pi{q^)[l + d^^\q^)] + ... = 1 + S^\q^) + d^^^q^) + ..., (5.4) 

where (i^^^(q^) = — (i(^^(<7^) — Pi(<7^)[l — Po(9^)]- Doing the third iteration, one further obtains 

d(q^) = 1 - P2{q^)[l + d^^\q^) + d^'^H^q^)] + ... = 1 + d^'H?') + ^^'H?') + d^^\l^) + -, (5-5) 

where d^^\q^) = —d'^^\q^) — d'^'^\q^) - P2{q^)[^ - Pi{q^){l - Pf){q^))], and so on for the next iterations. Thus up to 
the third iteration, one finaUy obtains 

°° A2 r A2 A2 t 

d{q^) = J2 d^"'\9') = 1 - i^2{q') + ^] 1 - [UAq') + ^][1 - Mq') ^] + - ■ (5.6) 

We restrict ourselves to the third iterated term, since this already allows to show explicitly some general features of 
such kind of the nonlinear iteration procedure. 

A. Splitting/shifting procedure 

Doing some tedious algebra, the previous expression can be rewritten as follows: 

d{q^) = [l-U2{q^)+n,{q')Il2{q^)-Uo{q')Il^{q^)Il2{q^) + ...] 

+ ^[Il2iq')Al + Ili{q')Al - no{q')Il,{q')Al - Ilo{q')n2{q^)Al - Il,{q')n2{q^)Al + ...] 

'i 

- ^[Iloiq')AlAl + n,{q')AlAl + Il2{q^)AlAl + ...] 

- -[^2-^2^+ ° 4 +-]' (5-7) 
qz qz q'i 

so that this formal expansion contains three different types of terms. The first type are the terms which contain only 
different combinations of I\-rn(q^) (they arc not multiplied by inverse powers of q^); the third type of terms contains 
only different combinations of (A^/g^). The second type of terms contains the so-called mixed terms, containing the 
first and third types of terms in different combinations. The two last types of terms are multiplied by the corresponding 
powers of Ijq^ . Evidently, such structure of terms will be present in each iteration term for the full gluon form factor. 
However, any of the mixed terms can be split exactly into the first and third types of terms by keeping the necessary 
number of terms in the Taylor expansions in powers of q^ for n„i((7^), which are regular functions at small q^ . Thus 
the IR structure of the full gluon form factor (which just is our primary goal to establish) is determined not only by 
the third type of terms. It gains contributions from the mixed terms as well. 
Let us present the above-mentioned Taylor expansions as follows: 

n™(g2) = n,„(o) + {q^/^?)Ii^^\Q) + {q^ / ^?fIl^X^) + a„(g'), (5.8) 

since for the third iteration we need to use the Taylor expansions up to this order (here /i^ is some fixed mass squared 
(not to be mixed up with the tensor index)). For example, the mixed term {llq^)Yl2{q^)A'{ should be split as 

^ii2{q^) = 4 [n2(o) + (//f,^)n^^Ho) + o(/)l = ^U2{0) + ain«(o) + o{q^). (5.9) 

Here and everywhere below am — (A'^/fi^), tti = 0, 1, 2, 3, ... arc the dimensionlcss constants. The first term now is to 

be shifted to the third type of terms and combined with the term (— 1/<7^)A2, while the second term ain2 {0) + O{q-^) 
is to be shifted to the first type of terms. All other mixed terms of similar structure should be treated absolutely in 
the same way. For the mixed term {—l/q*}Ilo(q'^)AlA2, one has 
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A?Ai 



g" 



.2/,,2^tt(1) 



no(0) + (g7M')n^^^(o) + [q' / n'Yn^,' {Q) + 0(9^) 



2/,,2a2tt(2)/ 



^no(o) ~ ^a2n«(o) - aia^n^^Ho) - 0{q^) 



(5.10) 



Again the first and second terms should be shifted to the third type of terms and combined with terms containing 
there the same powers of 1/g^, while the last two terms should be shifted to the first type of terms. 
Similar to the Taylor expansion (5.8), one has 



n™(g2)n„(q2) = n„„(g2) ^ n„„(o) + (gV^^^^w (o) + (^2/^2)2^^(0) + Omniq^) 

Then, for example the mixed term {—1 / q^)Ilo{q^)I{i{q^)/^\ can be split as 



(5.11) 



2 TT_ /„2\ 



A2 r 

^2 



(1)^ 



no(g^)ni(g^) = --f[noi(o) + ((77M'Kr(o) + o(g") 

= -^noi(0)-a2n^V(0) + 0(92), 



(5.12) 



so again the first term should be shifted to the third type of terms and combined with the terms containing the 
corresponding powers of Ijq^ , while other terms are to be shifted to the first type of terms. 

Completing this exact splitting/shifting procedure in the expansion (5.7), one can in general represent it as follows: 



d{q') = (A7g^)Bi(A,a,e,5') + (A7'?')'i?2(A, a, ^, g^) + {^' /q'fB:,{\a,^,g') + h{q') + 



(5.13) 



where we used notations (5.2), since the coefficients of the above-used Taylor expansions depend in general on the 
same set of parameters: A, 0,^,(72. The invariant function fsiq^) is dimensionless and regular at small q^; otherwise 
it remains arbitrary. The generalization on the next iterations is almost obvious. Let us only note that in this case 
more terms in the corresponding Taylor expansions should be kept " alive" . 

B. The exact structure of the general iteration solution 

Substituting the generalization of the expansion (5.13) on all iterations and omitting the tedious algebra, the 
general iteration solution of the gluon SD equation (3.21) for the regularized full gluon propagator (3.15) can be 
exactly decomposed as the sum of the two principally different terms as follows: 



A 2 '"^ 

D,,{q;A')^Dl,^^{q;A')+D^J{q) = ^^^,(g) ^(A2/q2)fe ^ $, ,„(A, a, C, 5^ 

^^ ' fc=0 m=0 



+ i 



T^,(g)^^„(g2)+^L^,(g)]J_ 



(5.14) 



where the superscript "INP" stands for the intrinsically NP part of the full gluon propagator. We distinguish between 
the two terms in Eq. (5.14) by the character of the corresponding IR singularities and the explicit presence of the 
mass gap (see below). Let us emphasize that the general problem of convergence of the formally regularized series 
(5.14) is irrelevant here. Anyway, the problem how to remove all types of the UV divergences (overlapping [ij] and 
overall [3, 0, fi 0) S Q ) i^ ^ standard one. Our problem will be how to deal with severe IR singularities due to their 
novelty and genuine NP character. Fortunately, there already exists a well-elaborated mathematical formalism for 
this purpose, namely the distribution theory (DT) [l^ to which the DRM [ll| should be correctly implemented (see 
alsoRef. [3). 

The INP part of the full gluon propagator is characterized by the presence of severe power-type (or equivalently NP) 
IR singularities ((j'^)~^~'^. A; = 0, 1, 2, 3, .... So these IR singularities are defined as more singular than the power-type 
IR singularity of the free gluon propagator (g^)~^, which thus can be defined as the PT IR singularity. The INP 



11 

part depends only on the transversal degrees of freedom of gauge bosons. Though its coefficients $fe_m(A,a,^, g^) 
may explicitly depend on the gauge-fixing parameter ^, the structure of this expansion itself does not depend on 
it. It vanishes as the mass gap goes formally to zero, while the PT part survives. The INP part of the full gluon 
propagator in Eq. (5.14) is nothing but the corresponding Laurent expansion in integer powers of q^ accompanied 
by the corresponding powers of the mass gap squared and multiplied by the sum over the g^-indcpendent factors, 
the so-called residues $fe(A, a,^,.g^) — X]m=o *^'=.'"('^''-'^'?'5^)- ^^'^ svan over m indicates that an infinite number of 
iterations (all iterations) of the corresponding regularized skeleton loop integrals invokes each severe IR singularity 
labelled by k. It is worth emphasizing that now this Laurent expansion cannot be summed up into anything similar to 
the initial Eq. (3.16), since its residues at poles gain additional contributions due to the splitting/shifting procedure, 
i.e., they become arbitrary. However, this arbitrariness is not a problem, because severe IR singularities should be 
treated by the DRM correctly implemented into the DT. For this the dependence of the residues on their arguments 
is all that matters and not their concrete values. The PT part of the full gluon propagator, which has only the PT 
IR singularity, remains undetermined. This is the price we have paid to fix the functional dependence of its INP 
part. In Refs. [lO, Il3| we came to the same structure (5.14) but in a rather different way. Concluding, it is worth 
emphasizing that both terms in the general iteration solution (5.13) are valid in the whole energy/momentum range, 
i.e., they are not asymptotics. At the same time, we achieved the exact separation between the two terms responsible 
for the NP (dominating in the IR {q^ -^ 0)) and the nontrivial PT (dominating in the UV {q^ -^ oo)) dynamics in 
the true QCD vacuum. This separation is unique as well, since for severe (i.e., NP) IR singularities there exists a 
special regularization expansion within the DT, complemented by the DRM, while for the PT IR singularity it does 
not exist [lllli|. 

VI. CONCLUSIONS 

Our consideration at this stage is necessarily formal, since the mass gap A^ remains neither IR (within the INP 
solution) nor UV rcnormalized yet. At this stage it has been only regularized, i.e., A^ = A^(A, a, £,, g^). However, there 
is no doubt that it will survive both multiplicative renormalization (MR) programs (which include the corresponding 
removal of both A and a parameters). For some preliminary aspects of the IRMR program see Refs. [iflllal- Anyway, 
how to conduct the UVMIlprogram is not our problem (as mentioned above, it is a standard one, and its description 
can be found in Refs. [3, 0, H 0, B M, llJ])- ^^ i^ worth noting that the mass gap which appears in the gluon SD 
equation cannot be in principle the same one which appears in the INP part of the general iteration solution, though 
we have identified them, for simplicity. 

It is important to emphasize that a mass gap has not been introduced by hand. It is hidden in the skeleton 
loop integrals, contributing to the gluon self-energy, and dynamically generated mainly due to the NL interaction of 
massless gluon modes. No truncations/approximations and no special gauge choice are made for the above-mentioned 
regularized skeleton loop integrals. An appropriate subtraction scheme has been applied to make the existence of a 
mass gap perfectly clear. Within the general iteration solution the mass gap shows up when the gluon momentum 
goes to zero. The Lagrangian of QCD does not contain a mass gap, while it explicitly appears in the corresponding 
gluon SD equation. This once more underlines the importance of investigation of the SD system of equations and 
identities [3, |j, [^, P| for understanding the true structure of the QCD ground state. We have established the general 
structure of the regularized full gluon propagator (see Eqs. (3.15) and (3.16)), and the corresponding SD equation 
(3.12) in the presence of a mass gap. 

In order to realize a mass gap, we propose not to impose the transversality condition on the gluon self-energy (see 
Eq. (3.5)), while preserving the color gauge invariance condition (3.14) for the full gluon propagator. This proposal 
is justified by the NL and NP dynamics of QCD. Such a temporary violation of color gauge invariance/symmetry 
(TVCGI/S) is completely NP effect, since in the PT limit A^ = this effect vanishes. Let us emphasize that we 
would propose this even if there were no explicit violation of the transversality of the gluon self-energy by the constant 
skeleton tadpole term. In other words, whether this term is explicitly present or not, but just color confinement (the 
gluon is not a physical state) gives us a possibility not to impose the transversality condition on the gluon self-energy. 
The existence of this term is a hint that the above-mentioned transversality might be temporary violated. Since gluon 
is not a physical state because of color confinement as mentioned above, the TVCGI/S in QCD has no direct physical 
consequences. None of physical observables in QCD will be directly affected by this proposal. For their calculations 
from first principles in low-energy QCD we need the full gluon propagator, which trasversality has been sacrificed 
in order to realize a mass gap (despite their general role the ghosts cannot guarantee its transversality in this case). 
However, we have already formulated a general method how the transversality of the gluon propagator relevant for 
NP QCD is to be restored at the final stage. In accordance with our prescription [131 it becomes automatically 
transversal, free of the PT contributions ("contaminations"), and it regularly depends on the mass gap, so that it 
vanishes when the mass gap goes to zero. 
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On the general ground (no truncations/approximations and no special gauge have been made) we have established 
the existence at least of two different types of solutions for the full gluon propagator in the presence of a mass gap. 
The so-called general iteration solution (5.14) is always severely singular in the IR (g^ -^ 0), i.e., the gluons always 
remain massless, and this does not depend on the gauge choice (this behavior of the full gluon propagator in different 
approximations and gauges has been earlier obtained and investigated in many papers, see, for example Ref. jlQ] and 
references therein). The massive-type solution (4.4) leads to an effective gluon mass, which explicitly depends on the 
gauge-fixing parameter, and it cannot be directly identified with the mass gap. Moreover, we were unable to make 
an effective gluon mass a gauge-invariant as a result of the renormalization, and therefore to assign to it a physical 
meaning. This solution becomes smooth at q^ ^ in the Landau gauge ^ = only. Both types of solutions are 
independent from each other and should be considered on equal footing, since the gluon SD equation is highly NL 
system. For such kind of systems the number of solutions is not fixed a priori. The UV behavior {q^ — > oo) of all 
solutions should be fixed by AF [3]. Due to unsolved yet confinement problem, the IR behavior {q^ — > 0) is not fixed. 
Only solution of the color confinement problem will decide which type of formal solutions really takes place. At the 
present state of arts none of them can be excluded. 

In summary, the behavior of QCD at large distances is governed by a mass gap, possibly realized in accordance with 
our proposal. The dynamically generated mass gap is usually related to breakdown of some symmetry (for example, 
the dynamically generated quark mass is an evidence of chiral symmetry breakdown) . Here a mass gap is an evidence 
of the TVCGI/S. Thus there is no breakdown of U{1) gauge symmetry in QED, since the photon is a physical state. 
At the same time, a temporary breakdown of SU{5) color gauge symmetry in QCD is possible, since the gluon is not 
a physical state (color confinement). In the presence of a mass gap the coupling constant becomes play no role. This 
is also a direct evidence of the "dimensional transmutation", g'^ — > A^(A, a,^,g^) d, [l^, 13, which occurs whenever 
a massless theory acquires masses dynamically. It is a general feature of spontaneous symmetry breaking in field 
theories. The mass gap has to play a crucial role in the realization of the quantum-dynamical mechanism of color 
confinement '511. 
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